Directing transport by polarized radiation in presence of chaos and dissipation 
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We study numerically the dynamics of particles on the Galton board of semi-disk scatters in 
presence of monochromatic radiation and dissipation. It is shown that under certain conditions 
the radiation leads to appearance of directed transport linked to an underlining strange attractor. 
The direction of transport can be efficiently changed by radiation polarization. The experimental 
realization of this effect in asymmetric antidot superlattices is discussed. 
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A great challenge for the future technology on micro- 
scopic scale is the ability to control transport created 
by external energy sources in presence of dissipation and 
noise. These sources may drive a system out of equi- 
librium and generate a transport which direction is re- 
lated to the system configuration in a rather nontrivial 
way. Such a directed transport appears in Brownian mo- 
tors, or ratchets, which now attract a great interest of 
the community [3, [3, At present the technological 
progress allowed to observe ratchets experimentally in a 
variety of systems including semiconductor heterostruc- 
tures vortices in Josephson junction arrays 0,01, cold 
atoms 7j, macroporous silicon membranes a\ and other 
systems. Biological applications of Brownian motors are 
also of primary importance Q, HI ■ 

The theoretical studies of ratchets have led to a number 
of interesting results (see 0] and Refs. therein) including 
such unexpected phenomena as absolute negative mo- 
bility when the transport is directed against an applied 
force 9]. However, up to now the studies of ratchets have 
been done mainly [iflj in the limit of strong dissipation 
neglecting the effects of second time derivative [lll^ or, in 
opposite, only the Hamiltonian dynamics has be analyzed 
Il2l | . In this Letter we consider an intermediate case 
with a moderate or weak dissipation in a two-dimensional 
dynamical system driven by monochromatic polarized ra- 
diation. Without dissipation the dynamics is completely 
chaotic and is characterized by a homogeneous diffusion. 
Surprisingly, the dissipation related to friction leads to 
appearance of a directed transport which direction can 
be changed efficiently by polarization of radiation. 

Our dynamical model describes the motion of noninter- 
acting particles on the modified Galton board in presence 
of friction and polarized monochromatic radiation. The 
board of rigid disks on a triangular lattice had been intro- 
duced by Galton in 1889 to analyze the dynamics of par- 
ticles elastically colliding with disks 0|. In absence of 
friction and radiation the dynamics is cornpletely chaotic 
as it had been proven by Sinai (see e.g. llj|). To break 
the symmetry of the Galton board we replace disks by 
semi-disks oriented in the x direction as shown in Fig. 1. 




FIG. 1: (color online) A chaotic trajectory (red/grey curve) 
on the Galton board of semi-disks (blue/black) at radiation 
with f — 4, 9 — tt/8, uj — 1.5 and dissipation 7 — 0.1. The 
trajectory starts near x — y = 0, the region shown corre- 
sponds to —27 < a; < —17, —10 < y < —5. Same trajectory 
on a large scale is shown in Fig. 3. 



Such a system has only translation and reflection y — » — y 
symmetries. 

We analyze the dynamics of particles on this modified 
Galton board under the influence of a linearly polarized 
radiation which gives a force directed by angle 9 to x- 
axis: f — /(cos 6*, sin 0) cos wt. In addition a particle ex- 
periences a friction force F/ = —7V directed against its 
velocity, and elastic collisions with semi-disks. A typical 
example of a chaotic trajectory in this regime is shown 
in Fig. 1. Without friction the dynamics is chaotic and 
a diffusive spreading in x and y directions takes place as 
well as a diffusive energy growth induced by ac-driving. 
However, due to the time reversibility, a directed trans- 
port is not possible in absence of friction. The introduc- 
tion of friction breaks the time inversion symmetry and 
creates a directed transport. 

To study numerically the dynamics of this model, we 
fix the disk radius r^ — \ and the particle mass m — 1 
so that the system is characterized only by the distance 
R between disk centers. In the following we consider the 
case of compact semi-disk board with i? = 2 (a moder- 
ate variation of R does not change qualitatively the dy- 
namics properties). Following the approach of H^, the 
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FIG. 2; (color online) Density distribution in the region 
—400 < X < —100, < 2/ < 600 averaged over 200 trajectories 
in a time interval <t < 2 - 10* for parameters of Fig. 1. Ini- 
tially all trajectories are distributed near x — y — with ran- 
dom velocity directions and — 1. For convenience the figure 
is clockwise rotated on n/2, density is proportional to color 
changing from zero (rose/black) to maximum (green/white). 

particle dynamics is simulated numerically by using the 
exact solution of Newton equations between collisions, 
and by determining the collision points with the rapidly 
converging Newton algorithm. This way at / = 0, 7 = 
the total energy is conserved with a relative precision of 
10~^^. A typical example of an average density distribu- 
tion obtained from a simulation of many trajectories is 
shown in Fig. 2. In spite of chaos the dissipation leads to 
an average directed transport of particles moving under 
angle ~ 2.3 to x-axis. 

The direction of transport depends on the polariza- 
tion of radiation as it is shown in Fig. 3 for typical 
trajectories computed on large time and space scales. 
The average velocity of directed flow can be written as 
V/ = vf{cosip,sintp). For radiation polarized along x 
[9 = 0) the average transport goes along —x direction 
(-0 — tt) while for radiation polarized along y [6 = 7r/2) 
the direction of transport is inversed {ij} = 0) ■ The fact 
that the particles move along x-axis for these 9 values is 
in agreement with the symmetry properties of the semi- 
disks board. However, in general the flow direction de- 
pends in a nontrivial way on the chaotic dynamics in 
presence of radiation and dissipation |^^. By changing 
the polarization angle 9 between and 7r/2 it is possible 
to change smoothly the angle of directed flow from tt to 
(see Figs. 3,4). The values oivf,ip and average velocity 
square v'^ — < + Vy > are determined from one long 
trajectory (up to times with ojt ~ 10^) or from averaging 
over 10 shorter orbits. The functional dependence 
found in this way is shown in Fig. 4. It varies moderately 
with system parameters but, in average, can be approx- 
imated by a linear relation -0 = vr — 29. We note that 
in presence of radiation the directionality properties of 
transport are much more rich compared to the case of 
static applied force where for disks the chaotic transport 
simply follows the force direction [islllij . 

To understand the properties of the directed transport 




FIG. 3; (color online) Directed transport for one trajectory at 
various polarizations of radiation 6 = 0, vr/S (same trajectory 
as in Fig. 1), 0.2l7r, 7r/4, 7r/2 (from left to right clockwise) at 
/ = 4,w = 1.5, 7 = 0.1. 

let us start from the analysis of its averaged character- 
istics. It is clear that in presence of chaos ac-driving 
leads to an effective heating increasing the particle ki- 
netic energy E — /2. This growth is stopped by en- 
ergy dissipation induced by friction. The balance of these 
two processes gives an average energy E of motion which 
due to chaos and dissipation takes place on a strange 
attractor This energy can be find on the basis of 

following estimates. Indeed, from the Newton equations 
of motion: E ^ vv ^ fvcosut. Hence, the diffusion 
rate in energy De = {AE)'^/t ~ E^t ~ pvl where the 
mean-free path Z ^ i? ~ 1 determines the collision time 
T ^l/v ^ RjE^I'^ and the space diffusion rate D ^ vl. 
Let us assume that the friction is weak and many col- 
lisions take place during the dissipative time scale I/7. 
Then, at the equilibrium the diffusion in energy is com- 
pensated by dissipation so that the average energy is de- 
termined in a self-consistent way as 

E = v'/2 ^ {DeH)'/' - {Ifhfl^ - {fhf" (1) 

where in the last relation we used that in our case 
I ^ R ^ 1. The same expression can be also obtained in 
a more formal way by writing the relation for velocities 
VmVn+i between two consecutive collisions n and 71 + 1 
from exact solution of the Newton equations. Then, as- 
suming that the dynamics is ergodic and chaotic on the 
energy surface, the squares of these velocities may be av- 
eraged over all directions and times between collisions. 
In the equilibrium both average and f^+i are equal 
that gives the relation The derivation of this ex- 

pression assumes that many collisions take place. This 
requires that the amplitude of oscillations induced by ra- 
diation a is larger then the distance between semi-disks: 
a = //(mw^) » R. Otherwise a trajectory drops on a 
simple attractor without scattering on semi-disks. 

The numerical data are presented in Fig. 5 and con- 
firm the dependence in almost 4 orders of magnitude 
range. The deviations appear only at large values of 7 
when the dissipation rate becomes comparable with the 
frequency of collisions. It is also interesting to note that 
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FIG. 4: (color online) Dependence of flow angle tp on polar- 
ization angle 6, measured in radians, for / = 4 and = 1.5 
at 7 = 0.2 (squares), 0.1 (circles), 0.05 (triangles) and 0.025 
(diamonds) and uj = 0.1 at 7 = 0.1 (x). Dashed curves are 
drown to adapt an eye. 



in agreement with ^ the dependence on radiation fre- 
quency UJ is rather weak. Indeed, a change of lu by an 
order of magnitude gives only a factor of 2 change in v'^ 
and in the flow velocity Vf (see Figs. 5,6). Even a larger 
change from uj = 2.5 to 0.025 at / = 4, 7 = 0.1 gives 
only 30% increase of and 50% drop of t;^ (data not 
shown). Such a weak dependence on lu can be attributed 
to finite time correlations, existing in the chaotic dynam- 
ics, which have been neglected in the derivation of 
The directionality of the flow is globally robust in respect 
to parameter variation (see Fig. 4). However, the depen- 
dence ip{0) has certain local changes for small 6 when ui 
becomes smaller or comparable with 7. 

The dependence of the velocity of the directed flow 
V/ is shown in Fig. 6 for 9 = 0. The numerical data are 
well fitted by the relation: 

Vf « ;2/3(^^/^2)i/3/^2 (7 > 7c) ; Vf « Ij/m (7 < 7^) 

(2) 

where 7c « (to///)^/^/40 and we assumed m = I = rd = 
1. The dependence of Vf on parameters for 7 > 7c re- 
minds the situation for a case of the Galton board with 
disks in presence of static force and friction (see 0,0|). 
Indeed, assuming that the flow moves in a given direc- 
tion its velocity can be estimated as Vf ^ fTc/m where 
Tc = l/v. According to 1^ this leads to (0) for 7 > 7c. 
However, in a difl'erence from 0, 0| the dependence is 
changed in the limit of small 7 <C 7c. Qualitatively, one 
may say that in this limit the particle velocity becomes 
very large, the chaotic collisions become very often and 
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FIG. 5: (color online) Dependence of average «^ on / and 
7 for 9 = at UJ = 1.5 and / =4(circles), 16 (squares), 64 
(triangles), 256 (diamonds) and 7 = 0.1 (*); the case at uj = 
0.1 and / = 4 is shown by (x); total interval of parameters 
variation is 5 • 10~^ <7<4, 2</< 256. The full line gives 
the dependence = (/^/7)^''"^/7. Logarithms are decimal 
here and in Fig. 6. 



this leads to a significant averaging and decrease of u/, 
which becomes proportional to the dissipation 7 that is 
responsible for appearance of directed fiow. With the 
variation of 9 from to 7r/2 the value oi Vf drops ap- 
proximately by a factor 2 {e.g. for data in Fig. 4). 

The results presented above show that the appearance 
of directed transport induced by polarized radiation is 
a robust phenomenon which is not very sensitive to the 
parameters of the model. Unfortunately we are not able 
to derive analytically the directionality dependence '4'i9). 
Qualitatively, we may argue that for 9 = the dissipa- 
tion leads to a trapping of particles in a funnel formed by 
semi-disks and the flow goes to the left. On the contrary, 
for 9 = Tr/2 the vertical oscillations induced by radiation 
push particles out of a funnel and in presence of dissi- 
pation the flow goes to the right. This indicates that 
approximately ijj = n — 29 in agreement with the numeri- 
cal data. However, it is desirable to have a more rigorous 
derivation in addition to these hand waving arguments. 
We also note that in the case of charged particles an in- 
troduction of a magnetic field leads to a reorientation of 
the directed flow following approximately the direction 
given by the average Lorentz force (we leave a discussion 
of this effect for further studies). 

The directed transport induced by polarized radia- 
tion demonstrates rather interesting and unusual prop- 
erties and it would be interesting to study it in experi- 
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The theoretical model discussed here should be also fur- 
ther developed to be more adapted for such experiments 
{e.g. Fermi surface effects and more realistic dissipation 
should be taken into account). 

We thank Kvon Ze Don for useful discussions. 
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FIG. 6: (color online) Dependence of the flow velocity Vf on 
friction 7 for 6 = and u> = 1.5 at / = 4, 16, 64 and 256 and 
Li; = 0.1 at / = 4 (same symbols as in Fig. 5). The straight 
lines show the asymptotic behavior Q at 7 > 7c (dashed) 
and 7 < 7c (full). 



ments with semiconductor heterostructures. The nowa- 
days technology allows to realize experimentally antidot 
superlattices (see e.g. 19]) and it is known that classi- 
cal chaotic dynamics plays an important role in trans- 
port properties of such systems (see e.g. H^]). How- 
ever, up to now only circular antidots have been stud- 
ied. In fact technologically it is rather simple to replace 
disks by semi-disks and to place the whole board in a 
microwave field. Moreover, the effects of microwave ra- 
diation on transport in disordered mesoscopic systems 
have been studied experimentally |2l| and it was shown 
that radiation may generate a directed transport in a 
case of broken space symmetry. Since disorder is sym- 
metric in average it is clear that the ratchet effect should 
be significantly larger in artificially prepared asymmetric 
structures as the Galton board of semi-disks discussed 
here. A typical set of experimental parameters can be 
rd ^ R ^ I ^ 1/im, Lj/iii = IGHz, m ~ 0.067me. For a 
microwave field of strength e 10 V/cm the amplitude 
of particle oscillations is a = ee/muP k, 60/im (a ^ R), 
so that the appearance of a simple attractor is excluded 
and a directed transport should appear due to presence 
of dissipative processes. A possibility to study transport 
in antidot lattices in presence of microwave radiation has 
been experimentally demonstrated in and therefor 
the experimental realization of our model should be pos- 
sible. Such experiments should allow to direct transport 
by radiation on mesoscopic scale and to study the dissi- 
pative processes in such systems from a new view point. 
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